Dynamics of quantum correlation for a central two-qubit coupled to an isotropic 

Lipkin-Meshkov-Glick bath 
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We investigate the behavior of quantum correlation, measured by quantum discord and entanglement, for 
two central spin qubits coupled to an isotropic Lipkin-Meshkov-Glick bath. We find that when the bath falls 
into different phases, both measurements of quantum correlation follow distinct paths. When the bath is in the 
symmetry-broken phase, as tuning the interaction between two-qubit and the bath being close to the interaction 
of bath spins, quantum discord and entanglement remain stable values with slight fluctuations. However, in 
the bath with the symmetric phase, quantum discord and entanglement always vary periodically with sharp 
fluctuations. Furthermore, the fluctuations can be enhanced by increasing the coupling strength between qubits 
and the bath. The critical point of quantum phase transition of the bath can be revealed clearly by the distinct 
behaviors of quantum correlation in the different phases of the bath. Besides, it is observed that quantum discord 
is significantly enhanced during the evolution while entanglement periodically vanishes. 

PACS numbers: 03.65.Ta, 03.65.Yz, ()3.67.Mn 



I. INTRODUCTION 



Nature utilizes correlation to bridge different things. En- 
tanglement, a kind of quantum correlation, has been taken to 
be fundamental resources for quantum information and com- 
munication 0], 01 ■ Extensive studies have shown that quan- 
tum entanglement is absolutely indispensable in such issues 
like superdense coding |[3|], quantum teleportation P], quan- 
tum cryptography Jst], etc. However, it was found that there 
are some types of nonclassical correlation other than entan- 
glement that also provide some advantages for the quantum 
systems over their classical counterparts Such a non- 

entanglement quantum correlation can be captured by quan- 
tum discord which is defined as the difference between 
quantum mutual information and classical correlation. Inves- 
tigations have shown that quantum discord is responsible for 
the computational efficiency of deterministic quantum com- 
putation with one pure qubit (DQCl) Isl-lioll. One feature 
of this intrinsically quantum correlation lies on that fact that 
quantum discord may have nonzero value for certain separa- 
ble states, while entanglement vanishes ||^[Tlt]. In general, 
for pure states quantum discord and entanglement entropy are 
equivalent, otherwise they are different. 

The difference between quantum discord and entanglement 
as figures of merit for characterizing quantum correlation is 
expected to be inherited and manifested in the phenomena 
of quantum phase transitions (QPTs) lfT2l - [l5ll . On one hand, 
since quantum discord grows towards the critical points in 
the context of QPTs, it work effectively as well as entangle- 
ment to signal QPTs. On the other hand, quantum discord 
may be more advantageous than entanglement at some con- 
ditions. For instance, pairwise quantum discord in infinite 
XY chain may be able to signal QPT for distances between 
sites where pairwise entanglement vanishes 1.14,1 . For XXZ 
spin chain at finite temperature, quantum discord spotlights 
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the critical point with respect to QPT while both entanglement 
and thermaldynamics quantity fails [15]. More recently, it was 
claimed that quantum discord is the best critical point estima- 
tor when one deals with finite-temperature systems fl^. 

For realistic quantum systems, the loss of coherence is in- 
evitable due to the interaction with the environment iITtIi . Re- 
cently, there have been increasing researches on quantum dis- 
cord dynamics under the effect of the environment [18-22J. 
In the Markovian environment, it has been shown that in all 
cases when entanglement suddenly disappears, quantum dis- 
cord vanishes only in the asymptotic limit [18]. Thus, it was 
argued that quantum algorithms based on quantum discord 
can be more robust than those based on entanglement. In the 
non-Markovian environment, the decline of quantum discord 
can be altered. [|20M22|l . 

Spin environment has been obtained much attentions ||23| - 
lls^]. Quantum coherence of a central two-qubit is rapidly de- 
stroyed as the environment being close to the critical point of 
its QPT ll26i I27I1 . Liu et al. have considered the effects of 
XY spin chain environment on quantum discord of a central 
two-qubit and showed that quantum discord can become min- 
imized close to the critical point of a QPT of XY chain BOll . 
Previously, Lipkin-Meshkov-Glick (LMG) model IsTll . origi- 
nally introduced in nuclear physics, has been proved to have 
a QPT 132-.34il . Purity of single central qubit coupled to an 
isotropic LMG bath witnesses the QPT critical point of the 
system Issll . Therefore it is intuitive to ask what would hap- 
pen for quantum correlation when a central two-qubit couples 
to such an isotropic LMG bath. 

In this paper, we will analyze the dynamics of quantum 
correlation of a central two-qubit, which interacts with an 
isotropic LMG bath. It is of interest to reveal the influence 
of the QPT on quantum correlation and the distinct behaviors 
of quantum coiTelation when bath is in two different phases. 
The paper is organized as follow. In Sec. [Ill we introduce 
our physical model and obtain the reduced density matrix of 
central two-qubit for two different phases. In Sec.|lIIl we cal- 
culate the quantum discord and entanglement, and present the 
main results. The last section is devoted to conclusions. 
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II. DYNAMICS OF A CENTRAL TWO-QUBIT COUPLED 
TO AN ISOTROPIC LMG BATH 



We consider two noninteracting central spin qubits coupled 
to an isotropic LMG bath. The total Hamiltonian can be ex- 
pressed as 



with 



Hm — 



(a ^ 0^ 
^ P K 
Q K y , 



(5) 



H — Hfi + Hi + Hs 



(1) 



with 



Hb = 

Hj = 
Hs = 



i<j '=1 

i=l 



a = (N^ -4M^)-2{M +1), 



J 



(2) 



2N 

~^N^_A{M + 2f)]-2{M+\), 

V 

-— ^2N{N + 2) - 8M(M + 1), 

~ V2A^(^ + 2) - 8(M + 1)(M + 2). (6) 



where He represents the Hamiltonian of the isotropic LMG 
bath. Hi denotes the interaction Hamiltonian between a cen- 
tral two-qubit and the bath, and Hs describes the central two- 
qubit. Here (7^,ia = x,y,z,i = 1,2, • ■ ■ ,A^) are the Pauli 
matrices of the i-th spin, cr"^^^ are the Pauli matrices of the 
central qubit a(b). N is the spin numbers of the bath, A is 
the coupling strength of the spins in the bath. A' the coupling 
strength between the central two-qubit and the bath. In the 
Dicke representation ll32l Issll . the above Hamiltonian can be 
rewritten as 



H = -^(StjS-^+S],S],-N)-2Sy^(S^S],+S-Stj)-2S„ 

(3) 

where 5^ = 5^ ± iS^ with S"^ = i = ^icr^ + 

cr-g) + j(cr^ + cTg) and Sr = 5(0-^ + cT^). The groundstate 
of isotropic LMG model lies in the subspace s pan ned by the 
Dicke states {|A?/2, M),M = -N/2,-- - ,A^/2) lHlll. It is 
easy to show that the energy corresponding to the Dicke states 
is A dependent, so the groundstate |G) is 1,34, .36.1 . 



|G> 



N N 
2 ' 2 



{0<A< 1), 
f,/(^)) (^>1), 



(4) 



where 1(A) denotes the integer nearest to N/{2A). It implies 
that there is a QPT at critical point A - 1. For the case 
of A > 1, the LMG bath is in a symmetric phase. For 
the case of < /I < 1, the LMG bath is in a symmetry- 
broken phase ll34l [36l l37ll . We below would study the dy- 
namics of quantum correlation for the bath in each phase. 
We simply denote triplets and singlet of two central spins by 

ii,i) = I mii,o> = -i=(i n> + i im ii,-i> = i u), and 



V2^ 

it)), respectively. 



10,0) = ^(1 Ti) 

In the invariant subspace Hm of H, spanned by the ordered 
basis {|A^/2,M)® |l,l),|A^/2,M-i- 1) » |1,0), |A^/2,M H- 2) ® 
|1, -1)), the total Hamiltonian can be expressed as a quasidi- 
agonal matrix with the blocks [i35l1 



The eigenvector corresponding to the eigenvalue Ej{j 
1, 2, 3) of Hm is directly obtained as 



I*?)) = PJ 



with 



M)®\l,l}+qj 



N 



M-t-l)®|l,0)-Hr; 



M+2)®\l,-l}, 
(7) 



Pj = 



1j = 





- a)HEj 


-y)^ + (Ej-y)^e 
(Ej - a)(Ej - y) 


+ (Ej 


— a)^K^ 




- a)\Ej 


-yf + (Ej-yy-e 
K(Ej - d) 


+ (Ej 


- af-K^ 




-a)HEj 


-yy + {Ej-y)^e 


+ (Ej 


- a^ie- 



(8) 



For simplicity, we consider the cases of Hm with three dis- 
tinct eigenvalues. Thus, (I*?))) is a set of orthonormal ba- 
sis of Hm- Then Hu is straightforward diagonalized as 
T^^HmT =diag{£'i, £2. £^3) with the matrix 



' P\ P2 P3 \ 

q\ qi qi 

n r2 r3 



(9) 



The dynamics evolution operator {/(f) - exp{-iHt) in the sub- 
space Hm can be obtained as llssll 



UM(t) = T 



,-iEit 



-iEjl 




PjqjC-'^'' 


-iEjt 






-iEjl 




qjrje-'^^' 



Uu Un Un 

U\2 U22 U23 
f/l3 f/23 U3J, , 



(10) 



It is worth noting that the Eq. ( fTOb is valid for the cases 
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-NI2 < M < N/2 - I, as {\N/2,M} ® \l,l},\N/2,M + 
1) i8> |1,0), \N/2,M + 2) » 1 1, -1)) is a three-dimensional in- 
variant subspace for these cases. When M = N/2 - 1, 
{\N/2,N/2 - 1) ® |1, l),|A^/2,A^/2> (g) |1,0)) forms a two- 
dimensional invariant subspace H^ii-i, and the dynamics evo- 
lution operator in this subspace can be calculated by the same 
procedure by setting a- = and y = in Eq. (|5]l. Besides, the 
states |A^/2,A^/2>® |1, 1) and \N/2,M) » |0,0> belong to one- 
dimensional invariant subspace, namely, they are the eigen- 
states of the total Hamiltonian H, the corresponding eigenen- 
ergies are -(A^-i-2) and 2AM^ -AN/2-2M, respectively. Thus 
the dynamics evolution of these states are exp[/(A^ + 2)t] and 
exp[-i(2AM^ - AN/2 - 2M)t]. 

To highlight the behavior of quantum correlation with re- 
spect to QPT, it is natural to assume the bath to be in the 
groundstate |G). Besides, we suppose that the central two- 
qubit is initially prepared in a class of X-structure state, 

Pab(O) = j{I + Etr^o-cTA ® '^P- '^^'"^ ^ identity op- 

erator on two-qubit, and (a = x, y, z) are real parameters 
to make pab{0) a legal quantum state. The evolution state of 
the total system is given by pioi(t) - U{t)p,ot{Q)U' {t), where 
Ptof(O) - \G){G\ ®pab{Q) is the initial state of the total system. 



where Uij,U'j, and U'j - 1,2,3) denote the matrices el- 
ements of the dynamics evolution operator in the invariant 
subspaces ///(^), and Hj^y2, respectively. 6(x) is the 

Heaviside step function, which equals to for x < and 1 for 
jc > 0. 



B. The LMG bath with symmetry-broken phase 

When the bath is in symmetry -broken phase (0 < A < I), 
the groundstate is |G) = \N/2,N/2}. As mentioned above, the 
dynamics evolution of the states |G) ® |1, 1) and |G) iSi |0, 0) 
can be directly obtained. |G) ® |1,-1> and |G) » |1,0) be- 
long to the invariant subspace //iv/2-2 and H^ii-x, respectively. 
One would note that the invariant subspace Hn/2-i is two- 
dimensional. By calculating the evolution state and tracing 
out the bath, we can obtain the reduced density matrix of the 
central two-qubit. 



Pabit) = Trbath[U(t)pUO)U'Ut)] 

with 



^(1 













(13) 



A. The LMG bath with symmetric phase 

When the bath is in symmetric phase (A > 1), 1(A) < N/2 
and the groundstate is |G) = \N/2, 1{A)}. One would notice 
that \N/2, M> ® 1 1 , 1 ), \N/2, M) » 1 1 , 0) and \N/2, M) » 1 1 , - 1 ) 
belong to three different invariant subspaces, Hm, Hm-\, and 
Hm-2, respectively. By the observation, we can apply Eq. (fTOl i 
to obtain the evolution state. The reduced density matrix of 
central two-qubit is obtained by tracing out the bath, and in 
the basis {| TT>, I ii>, I Ti>, I iT», it takes the form 



C 



1 + k- ~„ 1 -H fej- -H ky - k^ ~, ~ 
—r^il + \U'(,?) + -r^—\U'nt 



1 kjQ ky 



'221 



1 kx k\ kr 



Pabit) = TrhathWit^pam^if)] = 



with 



c 



J/ 



z 

z* c 



1 +k. 



4 

1 +k. 



-i\uur + \u[',n + 



if/i2r + if/2'3r) + 



1 kx ky k^ 



1 +L 



e(^!^-I(A)-ljun\^ + \U. 



33^ 



1 H" kx ~l~ ky . r\ 

1 +kr 



z = 



if/i2r + if/23r) + 

1 kx k^ 

8 ' 

IZJ^TT JT"* 

— - — UnU-,-, , 



S J/ 



(11) 



I + kjc + ky - k. 



\ + kx + ky - ky 



121 



22 1 



I + kjc + ky - ky 



U' 



221 



(12) 



~ _ kx ky 



i(N+2)tTTll* 

e ^33, 



(14) 



where U'-j{i,i - 1,2) and U"j{i,j = 1,2,3) denote the ma- 
trix elements of the dynamics evolution operator in invariant 
subspaces H^/i-i and H^ii-i- 

As shown above, all nonzero elements of reduced density 
matrices are time dependent. On contrast, when central two- 
qubit couples to an XY chain, the diagonal elements of re- 
duced density matrix remain constant |30]. This implies that 
effects of LMG bath on quantum correlation would be differ- 
ent from that of XY chain. 



III. QUANTUM CORRELATION OF A CENTRAL 
TWO-QUBIT 

Quantum discord is proposed to quantify all the nonclassi- 
cal correlation in the bipartite quantum system. For two-qubit 
state Pah, it is defined as difference between two versions of 
mutual information extended from classical to quantum sys- 
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tem 



QiPab) = KPab) - C(pab)- 



(15) 



Here I(pab) - S(pa) + S(pb) - Sipab) is the quantum mu- 
tual information which measures the total correlation in state 
Pah, where Sip) — -Tr(plogj p) is von Neumann entropy, 
and pa(h) is the reduced density matrix of pab- C(pab) as a 
quantifier of classical correlation is the maximum of quantum 
conditional entropy by performing one side measurement 



Cipat) = max ] 5 (pj - y ptS (p^'O 



(16) 



where {n['') denotes a set of von Neumann projectors on b, 
p^'* - T rti'i^l^ Pah^l^) I Pi is the state of a after obtaining out- 
come i on B, and p, = Trab(Tlf Pab^f)- 

Since the condition S{pa) - Sipb) is satisfied in expres- 
sion ( fTTI ) and ( fT3] l. the measurement of classical correla- 
tion assumes unique values, namely, C(pab) is irrespective of 
whether the measurement is performed on qubit a or b. We 
choose the complete set of von Neumann measurement to be 
of qubit b, i.e., {Ylf = = 1,2), with |1) = cos6»| T 

) + e'* sin e\ i) and |2) = e^'* sin 6\ T) - cos 0| J.). It is straight- 
forward to calculate the quantum discord of Eq. ( fTTT i [2Z1 . 



e(p„i) = min{Di,D2), 



(17) 



where 



and 

s o ^ 1-0, 1-0 1+0, 1+0 

D2^S(pa)-S (pab) — l0g2 — — logj , 

(19) 

with = Vc^^^^c)2T4ayiTizi)2. 

In order to compare the dynamics of quantum discord with 
the dynamics of entanglement, we choose the entanglement of 
formation as a quantifier of entanglement |38]. The entangle- 
ment of formation for two-qubit is a monotonically increasing 
function of concurrence |39], 




FIG. 1: (a) quantum discord and (b) entanglement as a function of 
the coupling strength A(A'^'') and the time t. Other parameters are set 
as k, = l,-ky = k, = l,andN = 500. 
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£ = -Alog2A-(l - A)log2(l - A), (20) 

where A = ^(1 + -y/l - with the concurrence ^. For the 
density matrix as Eq. ( fTTI ). one can directly calculate the con- 
currence as f = 2max{0, J/ - V^C, \Z\ - S]. By replacing 
the variables Ji, S, C, J/, and Z with ^, S, C, J/ and Z in 
expressions STfl and ( |20] | respectively, one can directly obtain 
the quantum discord and entanglement of formation for the 
density matrix in Eq. (fT3l l. 

Now, let us begin to analyze the dynamics of quantum cor- 
relation for the central two-qubit. By setting different values 
of ka{a = X, y, z), one can easily choose the initial state for the 



FIG. 2: Quantum discord (Q) and entanglement (E) as a function 
of the time / when the bath is in the symmetry-broken phase with 
A = 0.75, and (a) A' = 0.5, (b) A' = 2.0, (c) A' = 3.5, (d) A' = 5.0. 
Other parameters are set ask^; = 1, -ky = k^ = I, and A' = 1000. 



central two-qubit. In this paper, we will consider two cases, 
namely, the initial state to be pure and mixed. To investigate 
the dynamics behavior of quantum correlation, it is necessary 
to know the coupling strength between the cental two-qubit 
and the bath. First, let us study the case with A' - A. By set- 
ting kx = -ky - k^- 1, the central two-qubit is assumed to be 
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FIG. 3: Quantum discord (Q) and entanglement (E) as a function of 
the time t when the bath is in the symmetric phase with A = 1 .25, and 
(a) A' = 0.5, (b) A' = 2.0, (c) A' = 3.5, (d) A' = 5.0. Other parameters 
are set as kj, = 1, -fc,. = = 1, and N = 1000. 



FIG. 4: Quantum discord (Q) and entanglement (E) as a function 
of the time / when the bath is in the symmetry-broken phase with 
A = 0.75, and (a) A' = 0.5, (b) A' = 2.0, (c) A' = 3.5, (d) A' = 5.0. 
Other parameters are set as k^. = 1, -ky = k, = 0.2, and A' = 1000. 



in the Bell state -^(| TT) + I In Fig.[T] we plot the quan- 
tum discord and entanglement as a function of the coupling 
strength A and the time f. It is clear to see that the quantum 
discord behaves quite differently in each phase of the bath as 
A goes across A = I. Quantum discord remains as a constant 
unity in the symmetry-broken phase, while it oscillates with 
time t in the symmetric phase. The similar behavior of en- 
tanglement is also observed as depicted in Fig.lifb). Further- 
more, Fig.[T]also shows that entanglement periodically under- 
goes sudden death while quantum discord has local extrema 
at the same intervals of time. The common property of the 
oscillation of quantum correlation is that the frequency can be 
enhanced by increasing the coupling strength. The phenom- 
ena imply that the bath with symmetry-broken phase prefers 
to preserve the quantum correlation of the central two-qubit 
in the case A' = A. It is clear that the behavior of quantum 
correlation is quite different from that in XY chain ISOll when 
the environment approaches the critical point of QPT. 

To check the effects of the coupling strength A' on the quan- 
tum correlation, we do the numerical calculation for the ex- 
act expressions ( fTTj l and ( |20l l for fixed A values. Figure |2] 
shows the dynamical behavior of quantum discord and en- 
tanglement under different A' values when the bath is in the 
symmetry -broken phase with A = 0.75. As increasing the val- 
ues of A', the oscillation amplitudes of quantum correlation 
increase. In Fig. [3] we plot the dynamical behavior of quan- 
tum discord and entanglement when the bath is in the sym- 
metric phase with A - 1.25. The effects of the parameter A' 
becomes more significant in the symmetric phase than that in 
the symmetry-broken phase. The oscillation frequency of the 
quantum coiTelation apparently increases as increasing A'. In 
addition, it is clearly observed that quantum discord and en- 
tanglement synchronously oscillate. Fig.|2ja) also shows that 
the bath with symmetry broken phase preserves the quantum 
correlation for the case A' < A. Comparison between Fig|2]and 
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FIG. 5: Quantum discord (Q) and entanglement (E) as a function of 
the time t when the bath is in the symmetric phase with A = 1.25, and 
(a) A' = 0.5, (b) A' = 2.0, (c) A' = 3.5, (d) A' = 5.0. Other parameters 
are set as k, = 1, -yfc,, = k- = 0.2, and A' = 1000. 



Fig. [3] also shows that fluctuations of quantum correlation in 
the symmetry-broken phase are smaller than that in the sym- 
metric phase, which indicates the bath with symmetry-broken 
phase preserves a relatively stable quantum correlation. 

To study the dynamical properties of quantum discord when 
the central two-qubit is initially prepared in mixed state, we 
set the parameters - 1 , -^y - k, - 0.2, namely, the ini- 
tial state is Befl-diagonal state with pab(O) - 0.6|<1>+><<1>+| + 
0.4|>P+)<>I'+|, where 10+) = i(| TT) + I ii)) and \^'^) = i(| Ti 
) + I XT)) are Bell states. In Fig. |4] we plot quantum discord 
and entanglement as a function of the time t under different A' 
when the bath is in the symmetry -broken phase. As increasing 
the values of A' , the oscillation amplitude of quantum correla- 
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tion also changes with time t. Quantum discord and entangle- 
ment in this case oscillate in a more complicated way than that 
of the former cases illustrated in Fig.|2l Figure|5]shows the dy- 
namical behavior of quantum discord and entanglement when 
the bath is in the symmetric phase. The behavior of quantum 
correlation is quite different from the case of pure state with 
the same condition. In the cases of A' - 0.5 and 2.0, one can 
find that quantum discord gets maximum while entanglement 
gets minimum, and vice versa. As increasing the values of A' , 
the maximum values of quantum discord is also significantly 
enhanced while entanglement undergoes periodicity sudden 
death. This indicates the bath can induce quantum discord 
while destroyed entanglement at the same condition. 

IV. CONCLUSIONS 

In summary, we have investigated the dynamics of quantum 
discord and entanglement for a central two-qubit coupled to an 
isotropic Lipkin-Meshkov-Glick bath, which has symmetry- 
broken and symmetric phase. In the bath with symmetry- 
broken phase, quantum discord and entanglement of central 
two-qubit can remain stable values with slight fluctuations as 
tuning the coupling strength between the central two-qubit 
and the bath being close to the coupling strength of spins of 
the bath. However, in the bath with symmetric phase, quan- 
tum discord and entanglement vary periodically with sharp 
fluctuations. Furthermore, the fluctuations can be strength- 
ened by increasing the coupling strength between the cen- 



tral two-qubit and the bath, and the oscillation frequency of 
quantum correlation in the symmetric phase can be enhanced. 
Nevertheless, quantum correlation oscillate in the symmetry- 
broken phase with smaller fluctuations than that in the sym- 
metric phase. Our results imply that the bath with symmetry- 
broken phase prefers to preserve the quantum correlation in 
this sense. The critical point of quantum phase transition of 
the bath is clearly revealed by the distinct behaviors of quan- 
tum correlation in the two phases. In addition, in the case of 
central two-qubit with initial pure state, entanglement behaves 
as similarly as quantum discord except that entanglement oc- 
curs sudden death while quantum discord is still nonzero. In 
the case of central two-qubit with initial mixed state, quantum 
discord and entanglement become drastic. When entangle- 
ment suddenly disappears, quantum discord is significantly 
enhanced, which imply that quantum discord is more robust 
than entanglement under decoherence. Our results show that 
the isotropic LMG bath is able to provide an environment 
where the central two-qubit remains their coherence, and it 
may have potential application in quantum information pro- 
cesses. A further study on the effects of anisotropic LMG 
bath is expected. 
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